Stability of half-quantum vortices in + ipy superconductors 
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We consider the stability conditions for half-quantum vortices in quasi two dimensional px + ipy 
superconductor (such as Sr2Ru04 is believed to be). The predicted exotic nature of these excitations 
has recently attracted much attention but they have not been observed yet. We emphasize that 
an isolated half quantum vortex has a divergent energy cost in the bulk due to its unscreened spin 
current, which requires two half-quantum vortices with opposite spin winding to pair. We show 
that the stability of such a pair is enhanced when the ratio of spin superfluid density to superfluid 
density Psp/ps is small. We propose using various mesoscopic geometries to stabilize and observe 
these exotic excitations. 



The possibility of half-quantum vortices (1/2-qv's) in 
Px + iPy superconductors (SC's) has recently added a 
new 'spin' to the interest in such exotic SCs. The pre- 
diction that such vortices in (quasi-) two dimensional 
(2D) superfluids will have Majorana fermion zero modes 
bound at vortex cores, which render vortex statistics non- 
AbelianP,[3,i,|3, drives the excitement. A realization 
of excitations with such exotic statistics is of interest on 
its own right. Moreover, the possibility of exploiting non- 
abelian statistics for topological quantum computationfF] 
adds technological interest as well. 

1/2-qv's are topologically allowed in triplet superflu- 
ids due to their spin degrees of freedom. They were first 
sought after in thin films of^He-Afel: the best known ex- 
ample of Px-'ripy superfluid 7] . However, achieving a suf- 
ficiently thin parallel plate geometry is challenging and 
NMR experiments on "^He-A thin films failed to detect 
1/2-qv'sQ. Sr2Ru04 has recently emerged as a candi- 
date material for a quasi-2D (ab plane) PxMpy spin-triplet 
SC [13, U, [13], offering an alternate system to look 
for 1/2-qv's. Refs [3, [3] proposed experiments for probing 
and exploiting the exotic nature of the Majorana fermion 
core states of 1/2-qv's. Nevertheless, 1/2-qv's have not 
been observed in Sr2Ru04 [H, [3, Ell ■ For 1/2-qv's to be 
realized, it is crucial to carefully consider their energetics 
and schemes for stabilization. 

In this letter, we investigate the stability of the 1/2- 
qv's and discuss routes for their realization using meso- 
scopic samples. Vortex energetics in the context of 
Sr2Ru04 were first considered by Kee et al.^^ for 1/2- 
qv's in the he plane bound to a texture in the d vec- 
tor (Eq. ([1])). With the interest for the exotic statistics 
(which is only possible in 2D) in our minds, we focus on 
the ah plane vortices. Das Sarma et al. noted [^ that 
an out-of-plane field of the order of the spin-orbit decou- 
pling field[17j] can reduce the energy cost of 1/2-qv's in ah 
plane [3|. They further argued that the smaller vorticity 
of 1/2 qv's may energetically favor them over full qv's. 
However, as it was remarked in Ref.[3], an isolated 1/2- 
qv costs energy that is divergent in the system size due 
to its unscreenerf spin current. Indeed, fractional vortices 



that are allowed in multi-component SC generally tend 
to be energetically unfavorable [13] ■ Hence free isolated 
1/2-qv's cannot exist in the bulk. There are two ways 
to cut off this divergence: 1) two 1/2-qv's with opposite 
spin winding forming a pair in a bulk sample, 2) a 1/2-qv 
entering a mesoscopic sample. Through explicit energet- 
ics calculations, we find that 1/2-qv's can be stabilized 
when the ratio Psp/ Ps is sufficiently small. 

Stability Analysis - The order parameter of a triplet 



SC takes a matrix form in the spin space [12l . [1 



A(k) = 



A„(k) Atx(k) 
AiT(k) Au(k) 



id,, 



id,, 



The triplet pairing requires A 



(1) 

and the k depen- 
dent vector d(k) was introduced to parameterize the gap 
function. For each k, the unit vector d(k) represents 
the symmetry direction with respect to the rotation of 
Cooper pair spin. In a -I- ipy SC, the k dependence 
of d(k) is determined by the pairs having finite angular 
momentum projection m = I directed along 1. Hence 



d(k) = Aodexp(j(^j^j)]k x i| 



(2) 



where ip^^ is the azimuthal angle k makes around 1. For 
a (quasi-)2D system with 1 = ±z , Eq. ([2]) simplifies: 



d(k) = Aodexp(z(/?£), 



(3) 



for a single domain, where tp^ is the azimuthal angle k 
makes around the c-axis. Notice two independent contin- 
uous symmetries associated with d(k) : the spin rotation 
symmetry SO2 around d, and the U{1) symmetry com- 
bining a gauge transformation and an orbital rotation 
around 1 [20|. For fixed 1, we denote the 'orbital phase' 
for this combined U{1) symmetry by x- 

The spin degree of freedom represented by the d- vector 
allows for the existence of half-quantum vortices. In a 
singlet (single component) superconductor, 1/2-qv's are 
forbidden in order for the order parameter to be single 
valued. However in the Px + ipy SC, each component of 
the order parameter matrix can remain single valued by 
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simultaneously rotating the d-vector by tt while the or- 
bital phase X winds by tt. Since only % couples to the 
magnetic vector potential, this type of topological defect 
only carries a flux of hc/Ae, i.e. half a superconducting 
flux quantum $o = hc/2e. From Eq. ([1]), one can in- 
tuitively characterize a 1/2-qv as a vortex with a single 
unit of vorticity for one of the spin components and zero 
vorticity for the other [HI . 

In the absence of an external magnetic field, the spin- 
orbit (dipole) coupling favors alignment of d and 1, 
making 1/2-qv's energetically costly [23|- However, 
Knight-shift measurement on Sr2Ru04 by Murakawa et 
al. [17( suggests that a sufficiently large magnetic field 
H >200 G (< Hc2 750 G [13) along the c-axis may 
neutralize the dipolar interaction by fixing the d-vector 
orientation to the a5-plane (or Ru02 plane). While this 
opens the possibility of 1/2-qv's, the associated d-vector 
bending introduces 'hydrodynamic' spin terms in the gra- 
dient free energy. In the following, we will determine 
their effect on the stability of half-quantum vortices. 

The Ginzburg-Landau (GL) gradient free energy in its 
most general form in the London limit (i.e., the superfiuid 
density is taken to be constant outside vortex cores [20|): 

/grad = hpsVs - {Ps - ) (1 ' Vs)^] 
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[Psp(vJo'-(p.p-pL)a- vrf.)'] 



+ K"^"d,l,ndjln + aj{vsU\7 X Z)^. + — (VxA)^ (4) 
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where Vg is the superflow velocity, m is the fermion mass, 
and Ps, Ps and psp, pi'p are components of the superfiuid 
and spin fluid density matrix, respectively. Note that 
these are rank- two matrices, since there is one symme- 
try direction, 1, for the orbital degree of freedom. For 
a 2D single-domain, the fact that 1 is fixed leads to 
great simplification. The superflow velocity Vg is then 
{h/2m){\/x ~ 2eA/hc), as in a conventional s-wave SC. 
For vortex lines along the c-axis, Ps|| and Psp\\ are pro- 
jected out of the problem due to the absence of any vari- 
ation along the 1 direction (c-axis), and the superfiuid 
and spin current densities can be regarded as scalars. 
For d — (cos a, sin a, 0) in the a6-plane the free energy 
becomes 
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with the last term accounting for screening of the su- 
pcrcurrcnt, which is absent for the case of ^He-A thin 
film 6, 22]. 

Unlike a full qv, a single 1/2-qv costs a spin current 
energy that diverges logarithmically 



^sp 



(6) 
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FIG. 1; The energy difference (see the text) as a function 
of pair size, using ^ — 0.4A for Sr2Ru04 for different values 
of psp/ps- psp/ps = 0.3 (the blue solid curve), Psp/ps = 0.4 
(the purple dashed curve), and psp/ps =0.7 (the red dotted 
curve). The energy per unit length is in units of ^q/IGtv^X^. 



per unit length, where ^ is the core radius and R the 
lateral sample dimension. This divergence is due to the 
absence of screening, for the d-winding (the Va term in 
Eq. (O) does not couple to the electromagnetic field. 

A pair of 1/2-qv's with opposite sense windings in d 
would be free of such divergent energy cost. The energy 
per unit length of a pair of 1/2-qv's separated by ri2 is 
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where Kq is a modified Bessel function and X — -^^ 
the London penetration depth. The first term of Eq. ([7|) 
is the screened self-energy of superflow. The second term 
accounts for the magnetic interaction between the vor- 
tices, and the third term, in which the logarithmic diver- 
gence of Eq. ([6|) is canceled out by the interaction term, 
is purely due to spin flow. However, such a pair is topo- 
logically equivalent to a single full qv (see Fig. [TJa-b)) 
with the self-energy 



E 



full _ 



$2 
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Hence the most relevant question, regarding the stability 
of 1/2-qv's in the bulk of a SC is whether a pair of 1/2- 
qv's can be generated by the decay of a full qv. (Note 
that Eqs. ([7][5]) are approximately equal at ri9 =f .ll^ot 

We find that a 1/2-qv pair in bulk Sr2Ru04 will only be 
stable if Psp/ps satisfies certain conditions. In Fig. [l] we 
plot i;Jj^Jf(ri2)-£;f"" forri9> We used k = X/^^ 2.5 
based on Xab ~ 152 nm [l2,ll3l and the estimated value of 
the Ginzburg-Landau coherence length ^ab ~ 66 nm [l2| . 
For all three values of Psp/Ps shown, a finite equilibrium 
pair size requii > C exists. However, we note that while 
the London approximation predicts a minimum of the 
energy difference E^^l^{r 12) — i?^"" at a finite value ri2 = 
^'cquii for Psp/ Ps < 1, the London limit is only valid when 
Psp/ps is smaU enough to give requii > C 

We now discuss the ratio Psp/ps between the neutral 
spin superfiuid density Psp and the charged mass super- 
fiuid density ps in Eq. (O . In the case of '^He, this 
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ratio has been derived theoretically in terms of Landau 
parameters [v , 2^, 13, [l^ . Using sum rules for the longi- 
tudinal current-current correlation, Leggett showed that 
Psp < Ps [13 ■ This follows from the fact that the mass 
current is conserved even in the presence of interactions 
due to Galilean invariance, while the spin current is not, 
since a corresponding symmetry is absent. In the case 
of Sr2Ru04, the lattice breaks the Galilean invariance. 
However, since all interactions that scatter the mass cur- 
rent also must scatter spin current, but not vice versa, 
one still expects psp < Ps [23|. While the actual value 



of Psp/ps is unknown for Sr2Ru04, it has recently been 
measured to be ~ 0.3 near T = in ^He-A (2^. Since 
many Fermi liquid properties of Sr2Ru04 are similar to 
that of ^He Psp/Ps for Sr2Ru04 is also not antici- 
pated to be much larger than ~ 0.4. If so, the 1/2-qv 
pair would have a robust energy minimum (see Fig. [T]). 

Although the above analysis implies that a pair of half 
quantum vortices can be stable against combining into a 
single full qv, the vortex core energy can prevent a single 
full qv from decaying into a pair of half quantum vortices. 
Since the core interactions may be substantial for ri2 < 
it is not guaranteed that the finite equilibrium separation 
in our analysis corresponds to a global energy minimum. 
When a finite separation rcquii > configuration only 
represents a metastable state, its formation depends on 
the history of the sample. In this full vortex 

would not decay into two 1/2-qv's. If Psp/ps is so large 
that Toquii ^ ^, even a local minimum might be absent. 
However if Psp/ps ~ 0.4, as in the case of '^He, this is 
rather unlikely. 

Mesoscopic geometries - Given that the experimentally 
unknown core interaction may be si gnif icant, and neutron 
scattering ,13J and vortex imaging 3, 13] found no evi- 
dence for vortex dissociation in macroscopic samples, we 
propose utilizing mesoscopic geometries. A sample size 
of order A would promote the chance of a 1/2-qv enter- 
ing by cutting off the divergent spin current energy and 
reducing the effect of supercurrent screening. Here we 
discuss a thin slab and a cylinder as two examples. 

We find that a slab of thickness L — 2X (see Fig. [2]) 
will only permit the entry of 1/2-qv if Psp/ps < 
(2^/A)i4ri(2^/A). Consider the Gibbs free energy for a 
vortex at position x in the slab (see Fig. [2]), with wind- 
ing numbers rig and Ugp in x ^nd a, respectively, and 
a uniform applied field H: (rig = ngp = 1/2 and Ug = I, 
nsp = for a 1/2-qv and a full qv, respectively) 



G{x; ns,nsp) = s^nix; n^) 



+ 5]^(-ir-iebc(lx-a:^^l; n„n,p).(9) 

'5=TJ=l 

Here, Smix; Ug) is the energy (per unit sample thickness) 



due to the interaction with the Meissner current: 
e^{x; Us) H cosh(a;/A) - cosh[(L - 2^)/2A] 



'$o/47rA2 



cosh(L/2A) 

(10) 

The second line of Eq. ^ accounts for the boundary 
condition of vanishing current normal to each surface. 
We solve the boundary condition using an infinite set of 
image vortices at positions xj = (— 1)^ (x' =F jX) for j — 
l,...,oo with vorticities (— l)'''T.s,sp- Between a vortex 
with winding number {ns,nsp) in the slab and each of 
its image vortices at a distance r from it, there is an 
interaction energy (per unit sample thickness): 
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We plot G{x; 1,0) and G{x; 1/2, 1/2) in Fig.l^b). For 
a case when only the 1/2-qv is stable inside the sample. 
The energy difference between the two types of vortices 
at the center is estimated to be much larger than fcsTc 
at low enough temperatures (as large as W^ksTc/ fj.m at 
T ^ 100 mK). Hence, the entry of a 1/2-qv into the slab 
will be favored over that of a full qv. 

However, the fact that many vortices will enter the slab 
at different positions along its length complicates both 
further theoretical modeling and experiments. Thus, we 
also consider cylindrical samples, which can only accom- 
modate single vortices (of either type) if the radius is 
small enough [IJ]. For example, the fabrication of a 
sub-micron disk of Sr2Ru04, while challenging, is within 
reach of available technologies such as focused ion beam 
(FIB) milling. The entry of a vortex or fluxoid in a cylin- 
der or annulus geometry can be detected as an abrupt 
change in the magnetization as H is ramped up. As in 
the case of a thin slab, a 1/2-qv can be energetically fa- 
vored for small enough psp/ Ps- 

For simplicity, we only compare the free energies for 
fluxoids trapped in a thin, hollow cylinder. Although the 
fabrication of this geometry is more challenging, it has 
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FIG. 2: (a)A slab of thickness L = 2A = 0.304/im in an 
external field H. (b) G{x) (see Eq. (9)) in units of $g/(47rA)^ 
for a fuU-qv (blue solid curve) and for a 1/2-qv (red dashed 
curve) at H = 2V1G ~ 3.0$o/47rA^ for ^ = 0.4A and psp/ps = 
0.4. 
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an advantage over a filled cylinder or disk of eliminating 
the core. For a long, thin hollow cylinder of radius R and 
thickness d <C A, i?, the London approximation is nearly 
accurate [25j, and one can obtain the Gibbs free energy 
per unit length for fluxoids of either type in the pres- 
ence of an external field H from Eq. ([5|) via a Legendre 
transformation: 



spj 



^1/ 



= /3 



1 



Psp 2 
Ps 



$0 



(12) 

We defined (3 = dR/2X'^[26], and $ = ttR'^H, the exter- 
nally controlled flux through the cylinder. In the ground 
state, the values of Us and Ugp minimize G for a given 
value of <&. The n^p term in Eq. which is the 

only difference from the s-wave case, that allows for the 
1/2-qv possibility of Us = Usp = 1/2. At ^/^o = 1/2, 
G(1/2,1/2)<G(1,0) = G(0,0) if psp/Ps<(l+/3)"' and a 
1/2-qv will have lower energy than a full qv. The vortic- 
ity could be detected by monitoring the magnetization, 
as it was done for bila yer AI rings, which effectively form 
a two-component SC [27| . 

Conclusion - We considered the vortex energetics of 
Px + ipy SC in the London limit using a GL formalism. 
For an isolated 1/2-qv, we noted the importance of con- 
trolling a logarithmically divergent energy due to the un- 
screened spin current. While this divergence is regulated 
in a pair of 1/2-qv's of opposite spin winding, Psp/Ps 
needs to be sufficiently small in order to make a finite 
separation in such a pair (meta-)stable. Given the frag- 
ile nature of 1/2-qv's in bulk samples, we propose using 
mesoscopic geometries to look for 1/2-qv's. 

While for sufficiently small Psp/ps the pairs of 1/2- 
qv's can be energetically stable in bulk samples and a 
single 1/2-qv can be stable in mesoscopic samples, a 
number of real material features can affect the existence 
or complicate the observation of 1/2-qv's in Sr2Ru04. 
Such features include locking of the d-vector within the 
Ru02 plane, proliferation of microscopic domain walls, 
significant core energy effects, and boundary scattering. 
Nonetheless, given the exciting promise associated with 
these exotic excitations, our analysis forms a point of de- 
parture in the quest of realizing 1/2-qv's in Sr2Ru04 or 
other candidate materials for exotic superconductivity, 
such as UPts [111. 
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